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ABSTRACT 

The baryonic acoustic peak in the correlation function of galaxies and galaxy clusters pro- 
vides a standard ruler to probe the space-time geometry of the Universe, jointly constraining 
the angular diameter distance and the Hubble expansion rate. Moreover, non-linear effects 
can systematically shift the peak position, giving us the opportunity to exploit this clustering 
feature also as a dynamical probe. We investigate the possibility of detecting interactions in 
the dark sector through an accurate determination of the baryonic acoustic scale. Making use 
of the public halo catalogues extracted from the CoDECS simulations - the largest suite of N- 
body simulations of interacting dark energy models to date - we determine the position of the 
baryonic scale fitting a band-filtered correlation function, specifically designed to amplify the 
signal at the sound horizon. We analyze the shifts due to non-linear dynamics, redshift-space 
distortions and Gaussian redshift errors, in the range ^ z ^ 2. Since the coupling between 
dark energy and dark matter affects in a particular way the clustering properties of haloes and, 
specifically, the amplitude and location of the baryonic acoustic oscillations, the cosmic evolu- 
tion of the baryonic peak position might provide a direct way to discriminate interacting dark 
energy models from the standard ACDM framework. To maximize the efficiency of the bary- 
onic peak as a dynamic probe, the correlation function has to be measured in redshift-space, 
where the baryonic acoustic shift due to non-linearities is amplified. The typical redshift errors 
of spectroscopic galaxy surveys do not significantly impact these results. 
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1 INTRODUCTION 

Understanding the physical origin of the accelerated expansion of 
the Universe is one of the main drivers of the current cosmologi- 
cal research activity. The so-called A-Cold Dark Matter (ACDM) 
model is to date the most reliable cosmological scenario to explain 
such acceleration, due to its success in matching the observed prop- 
erties of the large scale structure (LSS) of the Universe. However, 
there are still some observational tensions at small and intermedi- 
ate scales, like the lack of luminous satellites in cold dark matter 
(CDM) haloes (see e.g. lNavarro et al.|[l996l : lBovlan-Kolchin et alj 
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Ettori 2003 ; McCarthy et 



ar yon tra 
al]|2007l) 



and the high velo cities detected 

in the large-scale bulk motion of galaxies (see e.g. IWatkins et alj 
l2009h or in systems of colliding galaxy clusters, as the famous 
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"Bullet Cluster" (see e.g. iLee & KomatsullioiOh . It is still unclear 
whether the above apparent disagreements are a consequence of 
some not well understood astrophysical phenomena at small scales, 
or if they really represent a signature of an underlying cosmological 
scenario different from the standard ACDM model. In any case, this 
provides an obvious motivation to explore alternative cosmological 
models that could explain these astrophysical anomalies. Further- 
more, the nature of the dark sector of the Universe, made of dark 
energy (DE) - the component responsible for the late time accel- 
eration - and dark matter (DM) - a kind of matter that negligibly 
interacts with radiation and can only be detected by its gravitational 
influence, is still unknown, even if it constitutes more than 95% of 
the total energy of the Universe. 

In the last fifteen yea rs, after the detection of the accelerated 
expansion of the Un iverse jRiess et al.ll998l : |Perlmutter et alj 19991: 
ISchmidt et al.lll99^ . widespread theoretical and observational ef- 
forts have been made to clarify whether the observed DE compo- 
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nent is just a cosmological constant A, or if its origin comes from 
other sources that can dynamically change in time. In this work, 
we investigate the LSS of the Universe in the so-called coupled 
(or interacting) DE models (cDE hereafter). This kind of cosmo- 
logical scenarios are based on the dynamical evolution of a clas- 
sical scalar field, (|), that plays the role of the DE and interacts 
with the CDM particles by exchanging energy-momentum (see e.g. 
I WetterichI 1 993: 1 Amendolall 200(1 12004 ) . Interestingly, cDE models 
can alleviate some observational tensions at small scales, that are 
presently not well und erstood in the ACDM cosmology, as men- 
tioned above (see e.g. iBaldi et aljboid ; iBaldi & Pettorindl20Ill : 
lLee&Baldill20Ilh . 

Baryonic Acoustic Oscillations (BAO) are the relic imprints 
left over by the primordial density perturbations in the early Uni- 
verse, when the photon pressure was counteracting the gravita- 
tional collapse of baryons, generating acoustic waves in the baryon- 
photon plasma. When the Universe became cool enough to allow 
the electrons and protons to recombine, at z ~ 1000, baryons and 
photons decoupled and the acoustic waves became frozen in the 
baryon distribution. Baryons then progressively fell into CDM po- 
tential wells and CDM was, in turn, al so attracted to baryon over- 
densities (see e.g. [Peebles & Yulll970l) . The BAO feature, corre- 
sponding to the scale of the comoving "sound horizon", can be de- 
tected in the matter distribution as a small excess in the number 
of galaxy pairs at ~100 /?^' Mpc , that creates either a single peak 
in the correlation function, ^(r), in configuration space, or a se- 
ries of peaks in the power spectrum, P{k), in Fourier space. As the 
sound horizon scale is well constrained by the cosmic microwave 
background (CMB) observations, the BAO location at different red- 
shifts provides a standard ruler to measure the expansion history of 
the Universe and to constra i n the DE equation of state, w(z) (see 
e.g. ISeo & Eisensteiijr2003l : lAlbrecht et aLll2006l : lAmendola et alj 
l20I2bl) . The BAO provide a powerful probe to investigate the na- 
ture of the late time acceleration, as they depend only slightly on the 
still uncertain astrophysical phenomena that could introduce sys- 
tematic bias in the cosmological constraints. 

Several independent measurements of the BAO feature in 
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planned to use these features for precision cosmol ogy, such as 
Euclicfl jLaureiis et alj|201 ll ; lAmendola et alj|2012hl) , BigBOS^ 
jSchlegel et al.l 201 ih . and the Wide-Field Infrared Survey Tele- 
scope (wfirstJI 

In order to exploit the BAO as a robust standard ruler, all the 
systematic effects have to be kept under control and the position of 
the BAO peak has to be measured with high precision. The tradi- 
tional estimators, ^(.?) and P{k), have some disadvantages that can 
affect the measure of the BAO at different redshifts. The former can 
be affected by the integral constraint at large-scale, if the cosmic 
number density of the population of mass tracers is not accurately 
known. The latter is instead affected by the uncertainty in the small- 
scale power, that is difficult to m odel due to non- linear effects. To 
reduce these possible svstematics. lXu et"aLl ( l201oh proposed a new 



robust estimator to analyze LSS and BAO, band-filtering the infor- 
mation in the two-point correlation function. The above estimator 
has been recently used to extract c osmological const raints from the 
BAO feature in the WiggleZ data l lBlake et alj|201 ll) . 

Non-linearities also affect the BAO as a standard ruler. Ac- 
cording to linear perturbation theory, the sound horizon scale 
imprinted in the early Universe remains unaltered as a func- 
tion of time. However, non-linearities cause a slight shift of 
the peak position relative to the linear predict i on (see e.g. 
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(RSD) further exacerbate this effect (see e.g. ISmith et alj|2008h . 
As a consequence, to calibrate the BAO as an unbiased standard 
ruler, the non-linear effects have to be accurately quantified in real- 
and redshift-space. Moreover, as we will discuss extensively in 
the next sections, the way that non-linearities and RSD affect the 
BAO feature depends on the assumed cosmological model. As a 
consequence, the BAO peak provides also a dynamical probe that 
could be of some help in discriminating between the ACDM model 
and alternative DE scenarios characterized by a different non-linear 
regime of structure evolution. 

In this paper, we will focus our investigation on configu- 
ration space and ext end to larger scales the work presented in 
iMarulli et alj ( l2012al) . analyzing the BAO fea ture in cDE mo dels 
with the band-filtered correlation function of IXu et aL 1 1 I2OIOI) . In 
particular, we focus on the shift of the BAO position due to non- 
linearities and RSD, in the redshift range ^ z ^ 2. Moreover, we 
quantify the impact of Gaussian redshift errors on the BAO fea- 
ture. To properly describe all the linear and non-linear effects at 
work, we employ a series of state-of-the-art N-body simulations 
for a wide range of different cDE scenarios - the COupled Dark 
Energy Cosmological Simulation project (CoDECS) jBaldi||2012cl) 
- that are used to investigate the impact that DE interactions can 
have on LSS, and in particular on BAO, with respect to the ACDM 
case. 

The structure of the paper is as follows. In <j2|and Sj3] we de- 
scribe the cDE models and the exploited set of N-body experiments 
analyzed in this work to simulate the LSS in these cosmologies. In 
5}4| we review our theoretical tools, describing the new statistic, 
(X)o(rs)> used in our investigation. Our results are described in SjS] 
where we analyze the CDM halo clustering at large scales, focus- 
ing on the shift of the BAO signal position through cosmic time. 
Finally, in Sj6|we draw our conclusions. 



2 THE COUPLED DARK ENERGY MODELS 

The cosmological scenarios investigated in the present work belong 
to the class of inhomogeneous DE models with species-dependent 
couplings (Damour et al. 19901). They have been indroduced as a 
possible alternative to the standard ACDM cosmology to explain 
the accelerated expansion of the Universe and to alleviate the "fine- 
tunin g" problems that affect the scenario with a cosmological con- 
stant lWetterich|[l993 ; lAmendolj[2000l) . These models are based 
on the dynamical evolution of a DE scalar field, (|), that interacts 
with the other components of the Universe by exchanging energy- 
momentum during cosmic evolution. In particular, the case of a 
coupling between DE and massive neutrinos has been proposed by 
lAmendola etai] ( l2008l) . while an interaction between DE and CDM 
particles has b een considered by e.g. IWetterichI lfl99^; lAmendolal 
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The BAO signal in cDE models 

Table 1. List of cosmological models considered in the CoDECS project and their specific parameters (see text for details). 
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l2012bh . Such DE-CDM interaction gives rise to a "fifth-force" be- 
tween CDM particles, mediated by the DE scalar field, that sig- 
nificantly modifies the gravitational instability processes through 
which cosmic structures develop, both in the linear and in the non- 
linear regimes (s ee e.g. [parrar & Rosen 2007; Ma ccio et alj|2004l ; 
iBaldi et al.l2oTol) . 

The background evolution of cDE cosmologies is described 
by the following set of dynamic equations: 



( lBrax&Martii]|l999h : 
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where an overdot represents a derivative with respect to the cosmic 
time t, H = d/a is the Hubble function, V{^) is the scalar field 
self-interaction potential, Mpj = l/VSltG is the reduced Planck 
mass, and the subscripts b,c\r indicate baryons, CDM and radi- 
ation, respectively. The "coupling function" 'r\c{<^) sets the strength 
of the interaction between the DE scalar field and CDM particles, 
while the quantity (t)r|£.((|)) indicates the direction of the energy- 
momentum flow between the two fields. If (i)r|(.((|)) > 0, the transfer 
of energy momentum occurs from CDM to DE, while the oppo- 
site happens when (t)'r|c((|)) < 0, due to the convention assumed in 
Eqs. m and This implies that the CDM particle masses, Mc, 
change in time, according to the following equation: 



dlnMc 
dt 



(6) 



Following lAmendoQ ( l2004t) . the coupling function can be ex 
pressed as: 



Tic 



:rioe ' 



(7) 



According to this equation, two different cases are possible: if 
Tji = 0, the coupling is constant, i.e. T|c((|)) = T|o; otherwise, the 
interaction strength changes in time due to the dynamical evolution 
of the DE scalar field (|) (see lBaldl201ll) . Both possibilities will be 
investigated in the present work. We will also consider two possible 
forms for the scalar field self-interaction potential, K ((])): 

• a "runaway" function given by the exponential potential 
dLucchin & MatarresdfTgsl ; IWettericljl 19881) : 



V{^) =Ae 



(8) 



• and a confining function given by the SUGRA potential 
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where the field (|) is expressed in units of the reduced Planck mass. 
The exponential model has a scalar field that always rolls down the 
potential with positive velocity (j) > 0, reaching the normalization 
value, (|) = 0, at the present time. The SUGRA model is instead 
assumed to start with the scalar field at rest in the potential min- 
imum (|)() = \fa in the early Universe, and the dynamics induced 
by the coupling allows for an inversion of the direction of motion 
of the field during the cosmic expansion history. This determines a 
"bounce" of the DE equation of state, w^, on the cosmological con- 
stant barrier, w,^ = — 1 , at relatively recent times, that significantly 
impacts the number density of massive CDM haloes at different 
cosmic epochs. Due to this peculiar feature, this kind of model has 
been termed the bouncing cDE scenario (see Isaldi 2012^ for a 
detailed discussion). 

The time evolution of linear density perturbations in the cou- 
pled CDM fluid and in the uncoupled baryonic component {^c.h = 
Spc./, IPc.b) is described at sub-horizon scales and in Fourier space 
by the following equations: 



-2H 



1-Tlf 



Sfc = -2/f 8/, + 47rG \pth + Pc^c 



5,+47iG[pi5/, + p,5,r,] , (10) 

(11) 



where, for simplicity, the dependence of the coupling function T|((|)) 
has been omitted. The factor Fc = 1 +4Ti^((j))/3 in Eq. l llOt rep- 
resents the fifth-force mediated by the DE scalar field, while the 
second term in the first square bracket at the right-hand-side of 
Eq. l llOt is the extra-friction on CDM fluctuation s arising as a con- 
sequence of momentum conservation (see e.g. lAmendolal l2004l : 



IPettorino & BaccigaluDiir2008l: iBaldi et alllioiol , 1201 ll) . A detailed 
description of the background and of the evolution of linear pertur- 
bations of all the cosmologies considered in this work can be found 
in ,Baldi. a012c.) and .Baldi. t2012ai) . 



3 THE N-BODY SIMULATIONS 

Similarly to all the other DE models proposed in the last decades, 
the cDE scenarios investigated in this work are barely distinguish- 
able from the ACDM model from their background expansion his- 
tory. Furthermore, despite the fact that such cosmologies determine 
a different growth of linear density perturbations as compared to 
the concordance ACDM scenario, such effect is highly degenerate 
with standard cosmological parameters like SIm and, most impor- 
tantly, CTg. Hence, to discriminate between these cosmologies and 
the ACDM one, it is useful to exploit also the non-linear regime of 
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Figure 1. The standai'd two-point correlation function of DM lialoes, ^(r) 
{upper panel), compared to the band-filtered coiTelation, (Bo(''.!) (lower 
panel). Both statistics are computed in the ACDM scenario, at z = 0. 
The green shaded bands mark the regions where the BAO feature is lo- 
calized: a peak at ~ llOft^' Mpcin the ^(i) correlation function, a dip at 
~ 130/!"' Mpc in a)o(r,). 



Structure fomation, which requires to rely on large numerical sim- 
ulations to accurately model all the non-linearities involved. The 
analysis described in the following sections is based on the pub- 
lic halo catalogues of the CoDECS simulations ( lBaldill2012cb . the 
largest N-body runs for cDE cosmologies to date. The range of 
cosmological models considered in the CoDECS project, which in- 
cludes the ACDM model as fiducial reference, and the constant, 
exponential and bouncing cDE cosmologies, are summarized in Ta- 
ble[T]along with their respective parameters. The viability of these 
models in terms of CMB observables has yet to be properly in- 
vestigated, in particular for what concerns the impact of variable- 
coupling and bouncing cDE models on the large-scale power of 
CMB anisotropies. Although such analysis might possibly lead to 
tighter bounds on the coupling and on the potential functions than 
the ones allowed in the present work, here we are mainly inter- 
ested in understanding the impact of cDE scenarios on the statis- 
tical properties of large-scale structures at late times, with a par- 
ticular focus on the role played by non-linear effects, and we de- 
liberately choose quite extreme values of the models parameters in 
order to maximize the effects under investigation. 

The simulation suite considered in this work is the L-CoDECS 
series, a set of collisionless N-body runs in a cosmological box of 1 
h^^ Gpcon a side, with 1024-' CDM particles and 1024-' baryonic 
particles, a mass resolution of mc{z = 0) = 5.84 x lO"^ Mq/A and 
nij, = 1.17 X lO'^ Mi7,//i for CDM and baryons, respectively, and 
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Figure 2. The survey geometry function <!>(;■) computed with Eq. )17t . The 
bin used for its computation is dr = 0.5 /i"' Mpc and Nr = SNo, where No 
and Nr are the number of haloes and random points, respectively. As can be 
seen, ^ i3?(r) ^ 1 and <3>(r) 1 when r^O. 

a force resolution of = 20 kpc/h. The simulations have been 
carried out with a modified version of the w i dely used parallel Tree- 
PM N-body code GADGET ( ISpringelll200^;lBaldi et al]|201(]h . 

The set of cosmological parameters at z = assumed in the 
CoDECS project are: Hq = 70.3kms"' Mpc"', £^cdm = 0.226, 
Q.UE = 0.729, Og = 0.809, Q.b = 0.0451 and % = 0.966, consistent 
with the "CMB -only Maximum Likehhood" results of WMAP7 
jKomatsu et alj|201lh . All the L-CoDECS simulations are normal- 
ized at z = zcMBi using the same initial linear power spectrum, 
and then rescaling the resulting displacements to the starting red- 
shift of the simulations, z; = 99, with the specific growth factor 
£)+ (z) obtained for each cosmological model by numerically solv- 
ing Eqs. {Tlinil- 

The CDM haloes have be en ident ified with a standard Friends- 
of-Friends (FoF) algorithm jDavis et al. 1985), using a linking 
length A. = 0.2 x d, where d is the mean interparticle separation. 
The results presented in this paper have been obtained using mass- 
selected sub-halo catalogues, composed by the gravitationally 
bound substructures that the algorithm SUBFIND dSpringel et alj 
I2OOII) identifies in each FoF halo. For all the simulations, we 
adopted the following mass range: M^i^ < M < M,^ax> where 
Mmin = 2.5 • IO'^Mq/A and M^ax = 3.6- 10■^l.l • 10'^4.9 • 
10i'',2.6- lO'^, 1.8- IO'^Mq/Zi at z = 0,0.55, 1, 1.6,2, respectively. 
In this paper, we only use mock subhalo catalogues and do not 
investigate which type of galaxies or clusters are more suited to 
achieve the highest accuracy and optimize our results. We defer a 
detailed discussion on this topic to a future paper. 



4 THE BAND-FILTERED CORRELATION FUNCTION 

The simplest way to analyze the spatial properties of the LSS of 
the Universe is to use two-points statistics. The probability, dP, of 
finding a pair with one object in the volume dVi and the other in 
the volume dV2, separated by a co-moving distance r, is given by: 

dP = ,i^[l+t,{r)]dVidV2, (12) 

where t,{r) is the two-p oint corre l ation function and n is the number 
density of the objects. IXu et al.l ( I2OIOI) introduced a new statistic, 
(£>l{rs), to analyze the LSS and, in particular, the BAO, that presents 
some good advantages with respect to the traditional estimators, 
t,{r), and its Fourier transform, P{k). In this paper, we will focus 
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Figure 3. The band-filtered correlation function (Bo(rj) for the different cosmological models considered in the CoDECS simulations, measured in real-space. 
The BAO signal is clearly seen as a single dip. Different panels refer to different redshifts, as labelled. 



on the monopole cO/=o(''j)! that is directly linked to the traditional 
t,(r) and dominates the signal with respect to all higher-order terms 
(see IXu et al.ll201ol for more details). The basic idea of this new 
statistic is to band-filter the acoustic information in the correlation 
function as follows: 



cOo(''j 



-4k 



r. dr 
Jo ts 



(13) 



The two-parameter filter function, W{r,rs), must be compact and 
compensated (f r^W{r,rs) dr = 0) with a characteristic scale rs 
jPadmanabhan et al]|2007h . A convenient choice is to use the fol- 
lowing analytical expression: 



lV(x)= 4x^(1 -x)(i-x) 0<x<l, 
W{x)=0 otherwise. 



(14) 



where x = {r/rsf' . Figure [T] compares the standard two-point 
correlation function of DM haloes (upper panel) with the band- 
filtered correlation {lower panel). Both statistics are computed in 
the ACDM model, at z = 0. While the BAO appears as a single 
peak, at ~ llO/i^' Mpc , in the traditional t,{r), it is instead a sin- 
gle dip, at ~ 130/i^' Mpc , in the new statistic cOo(''.s). It has been 
shown that coo(rs) combines the advantages of both the correla- 
tion function and the power spectr um approaches, at the same time 
jXu et alj|2O10l : lBlake et alj|201 ih . Indeed, (SX){rs) is insensitive to 
small-scale power, which is difficult to model due to non-linear 



effects, and it is also insensitive to large-scale power, where it is 
slightly affected by the integral constraint. It is also easy to mea- 
sure, as it only requires to weight the pair counts, without having 
to bin the data. If we adopt the natural estimator for the two-point 
correlation function: 



2 DDjnti) 

Ur,^) = ^rrr-. r-1, 



(15) 



RR{r,ij) 

where DD{r.ij) and RR{r.ij) are the number of galaxy-galaxy and 
random-random pairs with separation r, and line-of-sight angle 
arccos(/j), then from Eqs. l ll4t and il5\ we can convert the inte- 
gral given by Eq. l ll3t into the following expression: 



2 ^ Wjn/r,) 



(16) 



where Njj and «/> are the number and number density of objects, re- 
spectively, and <I>(r,/j) is a normalization function that encodes the 
geometry of the survey and depends on the number of data points. 
In general, ^ ^{r,ij) ^ 1 and when r./j — > 0, ^{r,p) 1. For 
simulations with periodic boundary conditions, ^{r,ij) = 1, as the 
volume is effectively infinite. The factor 2 appears because every 
pair is counted twice. The survey geometry function, <I>(r,ju), can 
be estimated by filling the survey volume with randomly distributed 
points, taking into account that their number has to be much larger 
than the number of data points, in order to keep the shot noise in the 
random pairs smaller than in the object pairs. Hence, the function 
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Figure 4. The band-filtered correlation function, coo(j'j), in the ACDM cos- 
mology, at z = and z ~ 2. Red stars show the measured (Bo(''j)' while the 
black lines are the best-fit polynomial model. The vertical dotted lines mark 
the position of the BAO scale. 
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Figure 5. The band-filtered correlation functions in real- (blue stars) and 
redshift-space (red stars), for the ACDM cosmological model at z = 1. 
Black lines are their fitting polynomials. The dotted (solid) vertical black 
line indicates the position of the BAO scale in real-space (redshift-space). 



<5(r,^) can be derived counting the random pairs, RR, and using the 
following equation: 

RR{r,ti) = 27tn£,A'or2*(r,^)drdA/, (17) 

where dr and d/j are the bins used to count the RR pairs in the 
survey. Notice that only the RR pairs, used to estimate the function 
<5(r,^), are measured in finite bins, while the data are not binned, 
unlike the traditional correlation function measurements. 



5 RESULTS 

5.1 Measuring the BAO scale in coo(f J ) 

To construct mock halo catalogues from the snapshots of the 
Codecs simulations, we locate virtual observers at z = and put 
each simulation box at a comovine distance corresponding to its 
output redshifts (see lMarulli et alj|2012bl : lBianchi et al.ll2012l for 
a detailed discussion). With the above method, each snapshot box 
represents a finite volume of the Universe, and the survey geom- 
etry function, <I>(r,/j), is not constant, as it would be if we mea- 
sured the distances relying on the periodic boundary conditions. 
However, our simple mock geometry guarantees that is 
almost independent of jA, as we have explicilty verified. We use 
Eq. l |17t to estimate this survey geometry function, generating ran- 
dom catalogues five times larger than the halo ones and counting 
the number of random pairs in bins of dr = 0.5 Mpc . The re- 
sult is shown in Fig. |2] We measure the band-filtered correlation 
function, coo(rj), using Eqs. l ll6t and l ll7t . and linearly interpolat- 
ing the function <5(r) that we computed as just discussed. Fig. |3] 
shows the real-space £00(^5) of DM haloes, as a function of red- 
shift (z = {0.00,0.55, 1.00, 1.61,2.01}), measured in the range 
10 < r5[/i"'Mpc] ^ 190. The black solid lines show the ACDM 
predictions, while the other curves refer to the cDE models of the 
Codecs simulations, as reported by the labels. At small scales, 
rj < 50 /7^'Mpc, we confirm the results found by iMarulli et alj 
( l2012al) using the standard two-point correlation function. More 
specifically, we find that, at z = 0, the small-scale halo spatial prop- 
erties predicted by the cDE models are very similar to the ACDM 
ones. Going to higher redshifts, it is easier to discriminate between 
our DE models that predict different shapes of C0o(r!)- Future ex- 
periments, such as Euclid, will be able to accurately constrain the 



spatial properties of high redshift objects, thus helping in discrimi- 
nating between alternative cosmological scenarios. 

However, in this paper we are more interested in the horizon- 
tal shift of the BAO peak due to cDE non-linear dynamics, as this is 
not degenerate with ag, nor with the linear galaxy bias. So we will 
focus on low redshifts, where this effect is larger. To accurately lo- 
calize the sound horizon scale in rj-space, with the band-filtered 
correlation function (i.e the scale corresponding to the minimum 
value of coo(rs))' we fitted the function (1)0(^5) around the dip us- 
ing a polynomial function. As we are only interested in the po- 
sition of the BAO feature, it is not necessary to model the shape 
of GiQirs) at all scales. Instead, we restrict our analyisis around the 
BAO dip: Itjo — r\< 5/?^' Mpc where r^Q is the scale corresponding 
to the minimum value of Q)o(rso)- In the above range, the shape of 
band-filtered correlation function can be accurately modeled with a 
polynomial function of third order: 

4 

Pi{>-s) = Y.ai.r[, (18) 

1=0 

where a/ are the coefficients of the polynomial. We estimate the 
values for all the cosmological models considered in the CoDECS 
simulations and at different redshifts, fitting coo(rs) with the model 
given by Eq. il8t and finding the local minimum of the best-fit 
polynomials. As a case study, Fig.|4]shows the band-filtered corre- 
lation function around the dip, measured in the ACDM halo mock 
catalogues at z = and z 2: red stars are the computed coo(rj), 
while the black lines show the best-fit polynomials. The vertical 
dotted lines mark the position of the BAO scale. As clearly shown 
in Fig. 151 the polynomial function is an excellent model for K)o(''i) 
around the minimum. 



5.2 The shift of the BAO scale 

5.2.1 The impact of non-linear dynamics 

The BAO position changes slightly through cosmic time, as can 
be directly observed using numerical simulations. This shift is 
mainly driven by physical effects, as non-linear structure forma- 
tion at large scales, that moves the position of the BAO signal to 
smaller scales relative to CMB-calibrated predictions. This occurs 
as a consequence of gravitational instability being non local (see 



© 0000 RAS, MNRAS 000, 000-000 



The BAO signal in cDE models 7 



Real-space 



ACDM 
EXP001 
EXP002 
EXP003 
EXP008e3 
SUGRA003 




1. 
1. 

„ 0. 

o 0. 
to 

^ -0. 

-1. 
-1. 

1. 
1. 

„ 0. 
o 0. 

« 

< -0. 

-1. 
-1. 



Real-space 




0.5 



ACDM 
EXP001 
EXP002 

\ 



EXP003 
EXP008e3 
•SUGRA003 



+ 



Redshift-space 




_L_ 



_L_ 



1.0 1.5 
Z 



2.0 



Figure 6. Evolution of the BAO scale in the band-filtered correlation func- 
tion, due to non-linearities for real- (upper panel) and redshift-space (lower 
panel) from z ~ 0.5 to z = 2. 



e.g. ICrocce & Scoccimarrj|2008l) . Therefore, at higher redshifts, 
where non-linearities have a small influence, the BAO peak in the 
two-point correlation function is located at larger scales than at 
lower redshifts. This effect can also be noticed in the band-filtered 
correlation function, with the difference that the BAO feature in r,- 
space, i.e r^g, is now located at smaller scales for higher redshifts, 
due to the effect of the filter, given by Eg. (114b. This can be clearly 
seen in Fig.|4l the position of the BAO scale changes in time, from 
fjo ~ 127.5 at z = to — 125 at z = 2. In r^-space, both the 
scale and the amplitude of the BAO feature increase as redshifts 
decrease. 



5.2.2 The impact of redshift-space distortions 

RSD induced by g alaxy peculiar mot ions further shift the position 
of the BAO scale jSmith et"al]|2008h . To construct redshift-space 
halo mock catalogues, we estimate the observed redshift of each 
object via the following equation: 



Zobs=Zc^ (1+Zc) 



(19) 



where Zc is the cosmological redshift due to the Hubble recession 
velocity at the comoving distance of the halo, V|| is the line-of-sight 
component of its centre of mass velocity and c is the speed of light. 
Fig.|5]shows the effect of RSD on the BAO feature, for the ACDM 
model at z = 1. Blue and red stars show the real- and redshift-space 
band-filtered correlation functions, respectively. The vertical dot- 
ted and solid black lines indicate the position of the BAO scale in 



Figure 7. The percentage difference between cDE and ACDM model pre- 
dictions of the evolution of the BAO scale in the band-filtered conelation 
function, due to non-linearities for real- and redshift-space from z ~ 0.5 to 
z = 2. 



real- and redshift-space, respectively. The amplitude enhancement 
of CO()(rj) in redshift-space is caused by the Kaiser effect and is di- 
rectly proportional to the linear growth rate of cosmic structures 
( lMarumetaDl2012ah . Here, we are more interested in the hori- 
zontal shift, as this is not degenerate with CTg, allowing to break 
the degeneracy between the effects arising as a consequence of the 
DE-CDM interactions and standard cosmological parameters. As 
shown in Fig. [5] the BAO scale is shifted by '-^ 2 /i^'Mpcfrom 
real- to redshift-space. 



5.3 The BAO scale as a dynamical probe 

The shift of the BAO scale due to non-linearities and RSD depends 
directly on the underlying cosmology. When the BAO feature is 
used as a geometric standard ruler, it is necessary to accurately 
correct for the shift, to avoid systematics when constraining cos- 
mological parameters. Most importantly, the correction should be 
different for each cosmological model to test. On the other hand, 
if the sound horizon scale is measured with sufficient accuracy, the 
BAO shift can provide a dynamical probe to discriminate between 
different cosmological models. 

The interaction between DE and CDM shifts the matter- 
radiation equality to higher redshifts, as CDM dilutes faster rela- 
tive toJhej\CDNl£as£_and there was more CDM in the past (see 
e.g. lAmendola et al. I l2012ah . Moreover, in cDE models the evolu- 
tion of linear and non-linear density perturbations is significantly 
altered, as compared to the standard ACDM scenario. Specifically, 
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Figure 8. Upper panel: The evolution of the difference between the sound 
horizon scale in redshift- and real-space for the ACDM scenario. Bottom 
panel: Percentage difference between cDE models and ACDM model pre- 
dictions obtained from the Eq. j20t . 
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Figure 9. The impact of redshift errors on the position of the BAO scale, for 
the ACDM model at different reshifts (z = 0.55 - solid line, z = 1 - dotted 
line, z = 2- dot-dashed line). 



as discussed in ij2] the DE coupling induces a time evolution of 
the mass of CDM particles as well as a modification of the growth 
of structures, determined by the combined effect of a long-range 
fifth-force, mediated by the DE scalar degree of freedom, and of 
an extra-friction, arising as a consequence of momentum conser- 
vation. At large scales, these effects change the location and the 
amplitude of the BAO scale. Therefore, the clustering properties of 
CDM haloes at the sound horizon scale, in particular the combina- 
tion of the BAO shifts due to non-linearities and RSD, might show 
specific signatures, allowing to distinguish a cDE Universe from 
the ACDM scenario. 



5.3. 1 The cDE BAO in real- and redshift-space 

Using our mock halo catalogues, we can quantify the effect of non- 
linear dynamics in the band-filtered correlation function, for the 
different cDE models of the CoDECS project. Fig.|6]shows the evo- 
lution of the position of the BAO scale, rjo, in real- (upper panel) 
and redshift-space (lower panel), as a function of redshift and cos- 
mological model. In both panels we can see that v^q increases, due 
to non-linearities, as redshift decreases. Comparing the upper and 
lower panel, we can see the impact of the RSD on the position of 
r^o for every cosmological model from z ~ to z '-^ 2, i.e. the shifts 
of rso to higher scales from real- to redshift-space as described in 

In Fig. [7] we show the percentage difference between cDE and 
ACDM predictions, Arso = 1 00 • (r^ccDE - rso,ACDM ) /fsO^ACDM , 
where Tso.acdm is the ACDM BAO scale, and Tso.cDE its counter- 
part for every cDE model. The percentage differences between the 
values of in ACDM and cDE models are quite small, < 2% in 
the whole redshift range. The EXP008e3 model predicts the largest 
difference, that rises from ~ 1% to ~ 2% at z ~ 1.5 and z ~ 2, re- 
spectively. The r^o difference in the SUGRA003 scenario is always 
positive (< 0.5%) and almost constant up to z ~ 0.5. 

The most interesting result of Fig.[6]is the amplification effect 
of RSD. While, in real-space, Arjo is small for all the cDE models 
considered (except for EXP008e3, at high redshifts), the predic- 
tions of the different cDE models are much more clearly distin- 
guishable in redshift-space. To better quantify this effect, in Fig. [8] 
we show the difference between the sound horizon scale in redshift- 
and real-space. RSD move the sound horizon at larger scales (as 



highlighted by the black line in the upper panel, that shows the 
ACDM prediction), an effect that monotonically increases at low 
redshifts. More interestingly, the effect of RSD is strongly model- 
dependent, as can be clearly seen in the lower panel of Fig. [8] that 
shows the percentage difference between Ar^o in ACDM and cDE 
models: 

»z irvr, ('^sO.cDE """sCcDe) ~ ('^O.ACDM "^sO.ACDm) ,-„, 

Ar,„ = 100 — — , (20) 

'-'sO,ACDM 'sO.ACDM^ 

where r^Q^j^pj^ and r^Q^^myj are the positions of the BAO fea- 
ture in redshift- and real-space respectively, in the ACDM scenario, 
while r^Q and r^Q are their respective counterparts for every 
cDE scenario. 

For constant coupling models, the RSD effect is directly pro- 
portional to the coupling strength, P: the larger is (3, the most the 
BAO scale is shifted to large scales. Indeed, the strongest effect is 
found for the EXP003 model, over the whole redshift range, for 
which the BAO shift due to RSD is up to ~ 90%. On the other 
hand, the impact of RSD is almost negligible for the EXPOOl and 
SUGRA003 models. We note that Arso can be negative for these 
models, but the effect is small considering the statistical uncertain- 
ties. 



5.3.2 The impact of redshift errors 

As we have seen in the previous sections, non-linear small-scale 
effects can "propagate" to larger scales, significantly shifting the 
sound horizon in a cosmology-dependent way. Hence, it is manda- 
tory to investigate how redshift errors can impact the BAO fea- 
ture. Here we consider only Gaussian redshift errors, typical of 
both spectroscopic and photometric galaxy surveys, and do not in- 
vestigate the impact of the so-called catastrophic errors (the ones 
caused by the misidentification of one or more spectral features). 
Indeed, pure catastrophic outliers with a flat distribution simply re- 
duce the amplitude of (£)Q{rs) at all scales, not affecting the position 
of fjQ. On the other hand, BAO shifts possibly induced by system- 
atic misidentifications of spectral features need to be estimated with 
mock galaxy catalogues, taking into account the specific observa- 
tional effects of the survey under investigation. 

We estimate the redshift of each mock halo with a modified 
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version of Eg. 1191 

Zobs=Zc+'^{l+Zc) + — , (21) 
c c 

where the new term Oy is the random error in the measured red- 
shift (expressed in km/s). To assess the impact of redshift errors 
on the position of r^o, we measure the band-filtered correlation in 
halo mock catalogues with the following Gaussian redshift errors: 
Ov = {0,200,400,600,800, 1200} km/s. In Fig|9l we show our re- 
sults only for the ACDM scenario, as the effect of redshift errors 
does not depend on the cosmological model. The impact of Gaus- 
sian errors results to be negligible for o,. < 600 km/s and, at high 
redshifts (z > 1.5), even for larger values of o,,. On the other hand, 
for z < 1 and > 600 km/s, the sound horizon systematically 
shifts to larger scales as redshift errors increase. Typical spectro- 
scopic galaxy surveys have quite accurate redshift measurements 
(C,, < 100 km/s), hence the systematic effect of redshift errors at 
the BAO scale is tiny in these cases. However, larger redshift errors, 
typical of photometric surveys, could introduce a non-negligible 
shift in the sound horizon scale, that should be accurately corrected 
when extracting cosmological constraints from the BAO feature. 
Interestingly, if we consider the required redshift error limits of the 
Euclid survey s, i.e. a- ^ 0.001(1 +z), which corresponds to 600 
km/s at z = 1 jLaureiis et alj|201 ih , the error on the location of rso 
is ~ 0.5%, i.e. small enough to distinguish cDE models from the 
ACDM scenario. 



6 CONCLUSIONS 

In this paper, we quantified the impact of DE interactions on the 
clustering of CDM haloes at large scales, using the band-filtered 
correlation function (£)o{rs). We used the public halo catalogues 
from the CoDECS simulations, a suite of large N-body runs for cDE 
models. We accurately estimated the sound horizon scale in our 
mock halo catalogues, fitting the function COo('".s) with third order 
polynomials around the BAO scale. We investigated the effects of 
non-linear dynamics both in real- and redshift-space, in the ranges 
10 < < 190 and ^ z ^ 2. The main results of our analysis are 
the following. 

• We detected BAO shifts in C0o(''j), relative to the linear pre- 
dictions, for ACDM and cDE cosmologies: non-linear effects shift 
the sound horizon to larger scales, an effect that monotonically in- 
creases moving to lower redshifts. 

• RSD significantly amplify the non-linear BAO shift, at all red- 
shifts. On the other hand, Gaussian redshift errors < 600 km/s have 
a negligible impact at the BAO scale. 

• The combined BAO shift due to non-linearities and RSD is 
strongly model dependent. Thus, if the BAO scale is accurately 
measured, it can provide a tool to distinguish between ACDM and 
cDE models. 

To conclude, besides the common use as a standard ruler, the 
BAO feature can be exploited also as a dynamical probe to disen- 
tangle between cosmological scenarios characterized by different 
non-linear dynamics. The BAO shift due to cDE dynamics is gen- 
erally small, but RSD significantly amplify the effect, expecially 
at low redshifts. Moreover, the shift of the BAO scale is not de- 
generate with ag and, consequently, with the other paratemeters 
degenerate in turn wi th c?s. hke the linear galaxy bias and the to- 
tal mass of neutrino s jLa Vacca et al .l2009l : lKristiansen et alj2"oiol ; 
iMarulli et alj|201 ih . Thus, the BAO feature provides a strong and 



unbiased probe for detecting new dark interactions, and can be of 
great help in removing the gg -degene racy of other in tependent ob- 
servables, like the h alo mass function jCui et al.l201^ and the clus- 
tering anisotropies jMarulli et alhoim 
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